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Abstract 

As gravity is a long-range force, one might a priori expect the Universe's 
global matter distribution to select a preferred rest frame for local gravita- 
tional physics. At the post-Newtonian approximation, two parameters suffice 
to describe the phenomenology of preferred-frame effects. One of them has 
already been very tightly constrained (|a2| < 4 x 1CP 7 , 90% C.L.), but the 
present bound on the other one is much weaker (|ai| < 5 x 1(T 4 , 90% C.L.). 
It is pointed out that the observation of particular orbits of artificial Earth 
satellites has the potential of improving the a\ limits by a couple of orders 
of magnitude, thanks to the appearance of small divisors which enhance the 
corresponding preferred-frame effects. There is a discrete set of inclinations 
which lead to arbitrarily small divisors, while, among zero-inclination (equa- 
torial) orbits, geostationary ones are near optimal. The main ai-induced 
effects are: (i) a complex secular evolution of the eccentricity vector of the 
orbit, describable as the vectorial sum of several independent rotations; and 
(ii) a yearly oscillation in the longitude of the satellite. 

PACS numbers: 04.80,+z, 95.40.+S, 11.30.Cp 
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I. INTRODUCTION 



The absence of preferred frames in local experiments (or local boost invariance) is verified 
everyday in high-energy experiments but is much more difficult to test for the gravitational 
interaction. Metrically coupled tensor-scalar theories of gravity (a la Jordan-Fierz-Brans- 
Dicke), including general relativity, do not predict any violation of this invariance (see e.g. 
[ ij and references therein). On the other hand, one expects the Universe's global matter 
distribution to select a preferred rest frame for the gravitational interaction if it is mediated 
in part by a long-range vector field (or a second tensor field besides the unique one postulated 
by Einstein) ||. Tests of the boost-invariance of gravity in localized systems are therefore of 
special interest to determine the field content of the gravitational interaction, i. e. specifically 
whether gravity contains, besides the standard tensor contribution and an often considered 
scalar one, extra contributions due to the exchange of a vector or a second tensor field. 

In the post-Newtonian limit, all preferred frame effects are phenomeno logically describ- 
able by only two parameters, ct\ and a 2 @] . [Note that the post-Newtonian approach assumes 
that all the fields contributing to gravity are massless, or at least have a range much larger 
than the size of the considered iV-body system.] These two preferred-frame parameters 
contribute non-boost-invariant terms in the Lagrangian, depending on the velocities v° of 
the bodies with respect to some gravitationally preferred rest frame. More precisely, the 
post-Newtonian Lagrangian describing the interaction between N spherical bodies reads 

L y = L/3,j,f] + L ai + L a2 , (1) 

L^ = ^-m^(l-K) 2 /c 2 ) 1/2 
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where the effective coupling constant G A b for the gravitational interaction of bodies A and 
B is 



G A b — G 

in which rj denotes the combination of parameters 
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r? = 4/3-7-3-ai + -a 2 , (6) 
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while E A ra,v = — ((7/2) f A f A d 3 xd 3 x'p('x)p('x')/\x. — x'| denotes the gravitational self-energy 
of body A. Apart from the contributions of «i and a 2 to Gab (via rj), Lp^ >ri is the usual 
boost-invariant Lagrangian involving the Eddington parameters j3 and 7, which are sufficient 
to parametrize the whole class of metrically-coupled tensor-multi-scalar theories ((3 = 7 = 
1 in general relativity, while these parameters can take arbitrary values in tensor-scalar 
theories). The parameters ot\ and a 2 contribute to the effects associated with a violation 
of the strong equivalence principle (Gab 7^ G) only through the combined parameter rj. 
Observational tests of the parameter rj are discussed in the literature (e.g. ||). In the 
following we concentrate on the effects associated to the velocity-dependent Lagrangian L ai 
ofEq. (|). 

It has been shown in Ref. @] that the close alignment of the Sun's spin axis with the 
solar system's planetary angular momentum yields an extremely tight bound on a 2 

\a 2 \ < 3.9 x lCT 7 (90 % C.L.) . (7) 

This limit on a 2 is much stronger than the existing limits on the other post-Newtonian 
parameters (3, 7 and a± ||. Present experimental bounds on the Eddington parameters are 
at the level @ 

|7-l|<3xl0~ 3 , \/3- 1| < 3 x 10~ 3 (90 % C.L.) , (8) 



while the limits on a± are about six times tighter []3],|6|-||. More precisely combined orbital 
data on the planetary system yield || 

«! = (2.1 ± 3.1) x 10" 4 (90 % C.L.) , (9) 

while the bound obtained in the strong-field regime by analyzing binary-pulsar data is com- 
parable || 

|e*i| < 5.0 x 10" 4 (90 % C.L.) . (10) 



Recent theoretical developments in tensor-scalar cosmological models |T0[ suggest that 
a natural level for \/3 — 1| and I7 — 1| is ~ 10~ 5 — 10~ 7 . The generalization of this result to 
classes of gravitational theories involving extra vector or tensor interactions has not been 
worked out, but, by analogy, one might expect the present limit on «i to be too weak to 
constrain at a significant level the participation to gravity of extra vector or tensor fields. 
It is therefore worth discussing experiments having the capability of improving the existing 
limits on «i down to the 10~ 5 — 10~ 6 level. 

The object of the present paper is to show that artificial Earth satellites offer very 
promising tools to improve the precision of measurement of a\ by, possibly, a couple of 
orders of magnitude. Indeed, we find that the appearance of small divisors can considerably 
enhance preferred-frame effects when the semi-major axis of the satellite's orbit and/or its 
inclination with respect to the Earth's equatorial plane take particular values. Section II 
is devoted to the secular evolution equations satisfied by the orbital elements, whose a\- 
dependent contributions derive from the disturbing function (0). As a first example of a 
small divisor which enhances preferred-frame effects, we consider in section III the simple 
case of equatorial orbits, and show that because of a competition between tidal forces and 
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the quadrupolar moment of the Earth, there exists an optimal value of the semi-major axis 
which maximizes the perturbations due to ot\. It turns out that geostationary satellites are 
nearly optimal. The case of non-equatorial orbits is studied in section IV. We point out 
that preferred-frame effects can be enhanced by arbitrarily large factors if the inclination 
is sufficiently close to one of six specific values, so that the accuracy of measurement of a\ 
becomes limited only by the finite duration of the experiment. 

When deriving bounds on a\ or a 2 [such as Eqs. (|7p,(]9D,(|I0D], it is necessary to make 
a definite assumption about the preferred rest frame entering the Lagrangians (|3]),(H). The 
standard assumption ||[||||| , that we shall take up in the present paper, is to choose the 
frame defined by the cosmic microwave background [this essentially means that the range of 
the putative extra vector or tensor field responsible for the violation of local boost invariance 
is at least of cosmological magnitude]. However, we shall see in section V.A that (somewhat 
weaker) bounds on a\ can be obtained without the need of such a hypothesis. Indeed, the 
'absolute' velocity of an artificial Earth satellite can be decomposed as v° at = v sat + v ffi + w, 
where v sat is the velocity of the satellite with respect to the Earth, v e the orbital velocity of 
the Earth around the Sun, and w the velocity of the Sun with respect to the gravitationally 
preferred frame. The ai-contribution (|3]) to the interaction term between the Earth and the 
satellite reads then 

at Gm e m sat 

L «i=-^ 2 (v© + w) ■ (v sat + v e + w) . (11) 

This involves a term proportional to v e • v sat which leads to observable effects, although they 
are 12 times smaller than those proportional to w • v sat , since %/c ~ 9.94 x 10~ 5 whereas 
w/c ~ 1.22 x 1CT 3 [|TTJ. Section V.B is devoted to the effects involving a.\V§ that are 
neglected in the rest of the paper. In particular, we show there that the term proportional 
to w-v e in Eq. (|II]) leads to sizable deviations of the angular position of the satellite, which 
could be used to constrain the value of a>i independently of the tests proposed in sections 
III and IV. We finally summarize our results and give our conclusions in section VI. 



II. SECULAR VARIATIONS OF THE ORBITAL ELEMENTS 

The preferred-frame effects described by the Lagrangians (|3]) and (§) are ~ aiv 2 /c 2 
smaller than the leading Newtonian interaction, and lead, over one orbital period, to very 
small deviations of the position of a satellite. Taking into account the present limits (|9T)- (|T0|) , 
one expects ai-induced displacements |5x| ~ «iw t> sat r esat /2c 2 ^ 75 /xm (for r esat ~ 2i? ffi ), 
too small to be observable with present or foreseeable techniques. Fortunately, some of 
these perturbations build up beyond one orbital period, and can thus be enhanced to an 
observable level if one waits for a sufficiently large number of periods. We restrict our 
attention to such effects in this paper, i.e. we concentrate on secular variations of the 
orbital elements. Concentrating on secular effects has also the advantage of freeing us from 
the coordinate ambiguities present in orbital-period effects fl2| . 



When taking into account only the Newtonian potential generated by a spherical Earth 
Gm©m sat /r esat , the satellite's orbit is determined by six constants of motion: its semi- 
major axis a, its eccentricity e, its inclination / with respect to the Earth's equatorial 
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plane, the angle f2 between a direction of reference (a, 5) = (0, 0) |13j and the ascending 
node, the angle u between this ascending node and the perigee, and finally the quantity 
a entering the mean anomaly £ = f Q n{t')dt' + er, where n = 2n/P = (Gm e /a 3 ) 1//2 is the 
orbital frequency. [Numerically P = (a/R^) 3 ^ 2 x 1.406 h, when using R & = 6.371 x 10 6 m, 
Gm 8 = 3.986 x 10 14 m 3 s -2 .] To help visualizing the meaning of the angle a, it is useful to 
note that, in the limiting case of a small eccentricity e, the mean anomaly I can be identified 
with the angular position of the satellite, so that a can be viewed as the angle between 
the perigee and the satellite at t — 0. Figure 1 summarizes this notation. It will also be 
useful in the following to define an orthonormal triad (a, b, c) , where a is directed towards 
the ascending node and c = a x b is in the direction of the orbital angular momentum, i. e. 
orthogonal to the orbital plane. (Note that this triad is not the one used in references |),14|, 
where a was directed towards the periastron.) 

We find it convenient to use the method of variation of the elements, as described for 
instance in Ref. [15], to derive the secular variations of a, e, J, a, u and Q. For more 
generality, let us consider a perturbed two-body system (m^mg), define M = uia + tub, 
Xa = tha/M, Xb = tub/M, and write the Lagrangian as L = Lq + MXaXbR where 
L contains, besides the kinetic terms, only the Newtonian potential between spherical 
bodies Gm^ms/r 'as, and R is the disturbing function containing all corrective terms due 
to asphericities, tidal forces, relativistic effects, . . . We shall not consider here the motion of 
the center of mass of the two-body system [Jj~6~| , but concentrate upon the equations satisfied 
by the elements of the relative orbit 
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The brackets in the right hand side of Eq. ( |12d| ) indicate that the a-differentiation is ef- 
fected keeping the mean anomaly fixed {i.e. ignoring the implicit a-dependence contained 
in J ndt). When working to first order in R, the orbital elements can be replaced by their 
constant zeroth order values on the right-hand-sides. Then the secular variations of the 
orbital elements are simply obtained by replacing in equations (0) R by its average value 
(R), computed over one unperturbed Newtonian orbit [More precisely, one takes the average 
over the angle £, keeping the other elements fixed]. It is straightforward to integrate the 
terms proportional to 1/c 2 in the two-body versions of the Lagrangians (0) and (|3|). De- 
composing the 'absolute' velocities as B = va,b + w cm, where w CM is the velocity of the 
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center-of-mass of the system (A, B) with respect to the preferred rest frame, we find the 
contribution 



ai ( GM\ 2 fw CM \ 2 



a 1 fGM\ / 2 [w CM /c,c,e 



GM\ 2 2 (2-i- (3 + 2 + a x X A X B 87 + 7 - - Aa x ) 



H^M (1-^ - T J- < 13 > 

where the square brackets denote the mixed product (wcm/c) • (c x e), the vector e being 
the Lagrange-Laplace (-Runge-Lenz) eccentricity vector with norm e, directed towards the 
periastron of body A. The first line comes from the term proportional to Wq M in L ai ([3|), 
and the second from the term involving wcm ' ( v a + v s)- The third line is responsible for 
the usual relativistic perigee advance and comes from the 1/c 2 terms of the Lagrangian (|l]), 
with extra contributions due to the • v B term of L ai . We have not included in Eq. (]T3|) 
the contribution of «2, relying on the limit (0) to consider that this type of preferred effects 
is already plausibly excluded. 

Up to here, we have considered a two-body system with arbitrary masses (m^ms) for 
generality's sake. In the particular case of artificial Earth satellites, one can however neglect 
m sat with respect to m e , and therefore write M ~ m e , Xa = X sat « 0, X B = X e « 1. 
Hence, we see from the last line in Eq. (|13D that the additional contribution to the relativistic 
perigee advance generated by a% is completely negligible. As for the first term on the right 
hand side of Eq. (pT3[), it only contributes to the evolution of the element a and will be 
studied in section V.B. In the next two sections we concentrate on the evolution of the other 
elements, generated by the second term on the r.h.s. of Eq. (|T3|). 

The contribution ( |13D should be added to the usual Newtonian multipolar and tidal per- 
turbations, notably the effect of the Earth's quadrupolar moment (J2) and the tidal forces 
due to the Moon (denoted as C) and the Sun (denoted as 0), which are the dominant 
ones. For analytical simplicity, we shall not consider any other Newtonian perturbation in 
the present paper. Our estimates of the measurability of a\ will correspondingly be only 
indicative. To get reliable estimates of the bounds on a\ which could be experimentally 
obtained, it would be necessary to perform detailed numerical simulations taking into ac- 
count all known Newtonian and general relativistic effects, as well as all modelizable sources 
of noise. [There exist in particular resonances between several multipolar moments of the 
Earth which have significant effects on long times.] 

The secular contribution due to the Earth's oblateness can be written down very easily 

Gm 9 R% 2-3sin 2 / 
= j2 ~^T- 4(1 - e 2 )3/ 2 ■ (14) 

The contribution of tidal forces is more complicated (even if one neglects the eccentricities 
of the Earth's and the Moon's orbits, as well as their inclinations with respect to the Earth's 
equatorial plane) because it involves explicitly the angular positions of the Moon and the 
Sun. However, it can be simplified very much if the variations of the orbital elements (|i2|) 
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are studied on longer time scales, more precisely if one averages over the orbital periods of 
the Moon around the Earth (~ one month) and of the Earth around the Sun (one year). 
The variations of the satellite's orbital elements after such averagings can then be derived 
from 



a 2 



((^Tides)) = N 2 — [2 + 3e 2 - 3 sin 2 / (1 - e 2 + 5e 2 sin 2 u)] , (15) 



where 



N 2 = ^ + ^- (16) 

r ©c r ©0 

[The notation N is reminiscent of the fact that n e = (Gm Q /r^ Q ) 1 / 2 is the orbital frequency 
of the Earth around the Sun; note however that (Gmc/r^) 1 ' 2 is not the orbital frequency 
of the Moon around the Earth, which would involve + mc instead of mcj Beware that 
the above expression fll5|) should not be used to derive variations with periods < one year. 
The secular variations of the satellite's orbital elements can now be easily derived from 



Eqs. (12). They tell us that the semi-major axis undergoes no secular variation ((d) = 
0), and that the change of the inclination is negligible in the limit of a small eccentricity 
[(dl/dt) = 0(aie) + 0(e 2 )]. We shall study the equations satisfied by e, ui and Q in sections 
III and IV, and see that they can be rewritten more compactly as a vectorial equation for 
the eccentricity vector e in the limit e<l. Section V.B will be devoted to the equation 
satisfied by a. 

Finally, it should be noted that in the main body of the paper we always deal with the 
motion of a satellite as seen in a (locally inertial) geocentric frame. The ai-dependent effects 
linked to the connection between such a geocentric frame and a global, barycentric one will 
be briefly discussed in the concluding section. 



III. COMPETITION BETWEEN TIDAL FORCES AND EARTH'S OBLATENESS 

FOR EQUATORIAL ORBITS 

The aim of this section is to show on the simple example of equatorial orbits how small 
divisors can enhance the preferred-frame effects associated to a\. [As said above, we neglect 
the eccentricities of the Earth's and the Moon's orbits, as well as their inclinations with 
respect to the Earth's equatorial plane.] The secular variations of the satellite's orbital 
elements are then given by equations fljjD , where (R) is the sum of Eqs. fll3|) , flip) and fllBD 
with 1 = 0. However, the angular position of the ascending node is no longer well-defined as 
/ — > 0, and equations ( |12e|) and ( |12i|) become formally singular. As is well known [15|, this 



singularity is fictitious and taken care of by considering the evolution of the angle ui = to + Vt 
between the direction of reference (a, 5) = (0, 0) and the perigee, which stays well-defined 
in the limit / — ► 0. One finds 

du\ (1-e 2 ) 1 / 2 d{R) 



u N - — — + O(aie) , (17) 



dt I na?e de 

ai n 2 a [w©/c, c, e 
T ~c e 
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where 



ujn = n 



Rn 



3 N 2 



(18) 



is the Newtonian perigee advance due to the Earth's quadrupolar moment and tidal forces, 
and w e = w + v e the velocity of the Earth with respect to the preferred frame, v ffi denoting 
as in section I the orbital velocity of the Earth around the Sun. [It is useful to quote here the 
numerical value of § n J 2 (R®/a) 2 = (R®/a) 7/2 x 2.02 x 10" 6 s" 1 « (R e /a) 7 / 2 x 2vr/(0.10 yr), 
which will appear again in the following sections as the characteristic frequency for the 
Newtonian variations of the orbital elements u and Q; see notably Eqs. ( |29| ) and (|31"|).1 
Equation (O) together with the one satisfied by e 



de 
dt 



'~4 



c 



c 



+ 0(a ie 2 ) 



(19) 



can be rewritten as a simple vectorial equation for the eccentricity vector e in the limit 

e < 1 



de 
~dt 



c x 



uo N e + (k + k @ ) x c + 0(«ie 2 ) + 0{e 2 ) , 



where 



-ax 



Aa 2 c 2 



w 



(20) 



(21) 



is a constant vector (beware that Ref. defines a vector k which is equal to twice this 
value) and 



Aa 2 c 2 



(22) 



a yearly varying one. In the present and following sections, we concentrate on the effects 
of the constant vector k, leaving to section V.A a study of the effects associated to K©(i). 
Since we are considering equatorial orbits in this section, the vector c orthogonal to the 
orbital plane is the constant unitary vector parallel to the Earth's polar axis, and therefore 
c x (k x c) is a constant vector, namely the projection of k onto the equatorial plane. Hence a 
ai-type violation of local boost invariance has the consequence of adding a constant forcing 
term in the time evolution of the eccentricity vector which, if it were alone, would secularly 
'polarize' the orbit in the direction of the equatorial projection of w. However, the Newtonian 
precession term ujn cuts off the build up of this polarizing term and deflects it by 90° in a 
gyroscope- like way. More precisely, the solution of Eq. ( p0|) can be written as the vectorial 
superposition 



e = e N (t) + e F 



(23) 



where ejv(i) is a vector of constant norm rotating with angular frequency un in the orbital 
plane (usual Newtonian perigee advance), and is the fixed polarizing contribution due to 
the preferred-frame effects we are studying, namely 
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c x k ,„ jN 
e F = ^— ■ (24) 

Note the factor (u5jv) _1 in Eq. fl2~4|) which is the first appearance of a small divisor. The 
solution (|23| ) is formally identical to the one found for binary pulsars in Ref. 0, the only 
difference being that ujn replaces the relativistic periastron advance. The same kind of 



polarizing terms in the time evolution of e has also been pointed out previously in [14], in 
the totally different context of equivalence principle violation in binary systems. 

In geometrical terms, Eq. (|23|) means that the eccentricity vector e(i) traces out, during 
its time evolution (after averaging over an orbital period), a circle centered around e^. If 
a large enough segment of this circle (say about a quarter) can be monitored during the 
experiment, it should be possible to measure the position of its center, i.e. e F , with about 



the same precision that individual measurements of the eccentricity vector [17] . In turn, the 
precision of the measurement of e is related to the precision a x with which one can measure 
the satellite's position: roughly one expects a e »s c x /a. It is convenient to work with 
quantities homogeneous to displacements. Therefore we shall measure the ai-perturbations 
of e in terms of 

_ ai_ |c x w| (Gm 9 ) 3/2 (r) r\ 
dp _ ae F - 2 j2GmffijR | a _ 5/2 + N 2 a5 /2 > W 

which represents (when e F < ex, as expected from e F < ct\ x 5.7 x 1CT 5 < 3 x 1CT 8 ) the 
amplitude of the secular change in the distance to the perigee p = a(l — |e|). From the 
arguments just given, equating 5p with the position measurement precision <r x should yield 
an estimate of the precision with which a% can be measured over a time span T/4, where 

2tt , . 

T = — ■ (26) 

UJ N 

The denominator of the expression (^5|) is a function of the semi-major axis a, and is 
minimized for the value 

^optimal = (2J 2 Gm (S Rl/N 2 ) 1 ^ 

w 7.66i? e w 4.88 x 10 7 m (27) 

(which corresponds to an orbital period ~ 30 hours). In Eq. (p7|), we have used the 
numerical data J 2 = 1.08263 x 10~ 3 , Gm^/R^c 2 = 6.97 x 10~ 10 , i? e = 6.371 x 10 6 m, and 
iV = 3.56 x 10~ 7 s _1 . The maximal value of bp (for equatorial orbits) is thus obtained for 
this value of the semi-major axis, and reads 

^ Gm§ C ^equatorial 

Pmax ~ ai R^ 6iV(2J 2 )V2 " 

w ot t x 2.54 x 10 5 cm , (28) 

where the norm w/c = 1.22 x 10~ 3 and the direction (a, 5) = (11.2 h, —7°) of w have been 
extracted from the results of COBE 
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Present technologies make it possible to measure the position of a satellite down to cr x ~ 
1 cm, by using either laser ranging (as for the Laser Geodynamical Satellites 'LAGEOS') or 
Global Positioning System (GPS) receivers onboard. The displacement (|28| ) could therefore 
lead in principle to limits on ot\ of the order of 4 x 10~ 6 , i.e. two orders of magnitude 
tighter than the present bounds (||) and (|T0|). It is to be noted that a geostationary satellite 
(a fa 6.62i? ffi , period ~ 23 h 56 min) is near optimal. It yields 5p ~ a± x 2.4 x 10 5 cm, 
leading us to expect a precision et\ ~ 4.2 x 10~ 6 . Clearly, numerical simulations taking into 
account all known Newtonian perturbations and modelizable non- conservative forces should 
be performed in order to get more realistic estimates of the precision which can be reached. 

This simple case of equatorial orbits already exhibits one of the characteristic behaviors 
of preferred-frame effects enhanced by small divisors: the duration of the experiment must 
be large enough to take advantage of the small divisor effect. For instance, the time span of 
the experiment should be (at least) T/4 = k/2lon ~ 15.4 yr in the case of the optimal orbit 
(f2"7|), and ~ 12.4 yr for a geostationary satellite. 

Let us conclude this section by a comparison of the result (^) with the one corresponding 
to the best drag-free satellite we know: the Moon itself. Still neglecting the inclination of its 
orbit and that of the ecliptic with respect to the Earth's equatorial plane, we can repeat all 
of the above discussion, except that in Eq. flTE| ) only the tidal forces due to the Sun should 
be taken into account. Actually, higher-order terms in these tidal forces are not small in 
the particular case of the Moon (because its orbital period is non-negligible as compared 
to one year), and the coefficient N 2 of Eq. (|l"5|) should be replaced by 2.043 instead 
of merely n\ = Gm Q /r^ Q to yield the correct magnitude of ui N |15]| . The amplitude of 
the secular oscillation of the Moon's perigee distance is then given by the same formula 
(p5|) as above with a = r eC = 60.3 R®. The Earth's oblateness gives then a negligible 
contribution to the denominator of 5p in Eq. (|25|) , and one finds Spc ~«ix 4.5 x 10 3 cm, 



with period T m 8.72 yr. This effect was first discussed (in a different guise) in Ref. ||T8 
Since Lunar Laser Ranging (LLR) allows one to determine the Earth-Moon distance within 
2 or 3 centimeters, we expect that an analysis of LLR data could at best limit ot\ at the 
~ 5 x 10~ 4 level, i.e. not better than the existing bounds (p|) , fllCf) . However, as this 
would represent a new, independent test of preferred-frame effects, it would be interesting 
to perform explicitly a multi-parameter fit of LLR data including the contributions of a± 
(and for completeness a 2 , /3, 7) to the Lunar equations of motion. 



IV. NON-EQUATORIAL ORBITS AND ARBITRARILY SMALL DIVISORS 

In this section, we shall see that non-equatorial orbits (of any altitude) can give rise 
to arbitrarily small divisors if the inclination is near some special values. For reasons that 
will appear clear below, we concentrate upon rather low satellites, for which tidal forces are 
negligible compared to the influence of the Earth's quadrupolar moment. For this reason 
we do not take into account the disturbing function (|T5| ) in the present section, in order 
to avoid unnecessarily technical calculations. However, we show at the end of this section 
how the equations satisfied by the orbital elements can be solved when tidal forces are not 
neglected. 

The main difference between equatorial and non-equatorial satellites is that the orbital 
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plane is no longer fixed when 1^0. Indeed, equation ( |12f|) shows that it is precessing with 
a constant angular frequency 



(O) = -1 n J 2 f^j cos / + O(aie) + 0(e 2 ) + 0(N 2 ) , (29) 

which is not modified by a cui-type violation of local Lorentz invariance in the limiting case 
of a small eccentricity. We shall drop the angular brackets in the following, and denote 
this constant precessing velocity simply by Cl. Like in the previous section, the equations 
( |12b| ) and ( |12ej ) satisfied by e and u can be rewritten as a simple vectorial equation for the 
evolution of the components of e with respect to the vectors (a, b), which are part of a frame 
[defined at the beginning of section II and in Fig. 1] which rotates around the Earth's polar 
axis with the angular velocity (p9|): 

(d'e/dt) = c x [uj N e + (k + #c ffi ) x c] + 0(«ie 2 ) + 0(e 2 ) + 0(eN 2 ) , (30) 

where the prime in d'/dt denotes a time derivative in the rotating frame (a, b, c) and where 

u N = ~ n J 2 f^j 2 (4 - 5 sin 2 I) + 0(e 2 ) + 0(N 2 ) (31) 

is the Newtonian perigee advance due to the Earth's quadrupolar moment, k and being 
the vectors defined in Eqs. (pID , (p2[) . As in section III, we shall neglect in this section K e 
with respect to k. Since the projection c x (k x c) of k onto the orbital plane is no longer a 
constant vector because of the precession Q, the solution (^) of the previous section is not 
valid for non-equatorial orbits. However, it is easy to solve equation (|30|) which is just an 
inhomogeneous linear differential equation in e (in the limiting case of a small eccentricity). 
Let us look for a solution of the type e = e^v + Yli e ii where ejy is the usual (constant-norm) 
Newtonian eccentricity vector rotating in the orbital plane (a, b) with angular frequency 
u>n, and where the e^'s are some constant-norm vectors rotating in the orbital plane with 
constant angular frequencies Ui to be determined. The time derivative of e then reads 



d'e/dt = 




This has precisely the form of equation (|30"D, with k x c having been decomposed as a sum 
of constant-norm vectors (uJi — uJN) e i rotating with constant angular frequencies uoi in the 
orbital plane. Since k x c has a periodic motion in the orbital plane (with period 2ir/Cl), it 
obviously admits such a decomposition: 

k x c = -ai ^§ w (K+ + K + K_) , (33) 

where 
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I i COS I 

K± = cos 5 [=Fsin(fi — a) a — cos(fi — a) h] , (34a) 

K = sin S sin I a , (34b) 

(a, 5) denoting the right-ascension and declination of w. K± are rotating with angular 
frequencies =|=fi in the orbital plane, whereas Ko is a constant vector in this plane (directed 
towards the ascending node). The solution of Eq. (|30|) can therefore be written simply as 

e = e N + e + + e + e 

Gm m ( K+ K K 



4a 2 c 2 \u N + n uj n uj N -VtJ 
This generalizes the solution (|23|) obtained above for equatorial orbits. [In the limit I — > 0, 



the term involving K + becomes equal to the vector ep of Eq. (|24|) , in which Cjn + Vt 
was denoted ujn-} In geometrical terms, Eq. fl35|) means that the eccentricity vector e(t) 
undergoes a kind of epicyclic motion in the (a, b) plane: it moves (with angular velocity ujn) 
along a circle of radius |ejv | whose center moves itself (with a non-uniform angular velocity 
determined by Q) on an ellipse (e + (t) + e_(t)) centered around e . Hence we can now 
distinguish three different contributions to the «i-induced secular oscillation of the distance 
to the perigee: 

Gm ffi r (l±cos/)/2 

op± = ae± = oli — wcoso = 

4ac 2 to N ±Vl 

(l±cosJ)/2 

= Acoso K , (36a) 

4-5 sin 2 / =f 2 cos / 

. Gm ffi . sin J 

opo = aeo = cxi 5- w sm 

4a<r ujn 

= A sin 5 S1 - n/ 2r , (36b) 
4 — 5 sm 1 



where 



^4 = 0;! 



a \ 5/2 ( Gm e \ 1/2 _Re w 
R e Uec 2 y 3J 2 c 



w ai (a/i? e ) 5/2 x 6316 cm. (37) 

Figure 2 displays the amplitude of these displacements as functions of the inclination. They 
clearly exhibit poles (or 'resonances') for six particular values of/: 5p + diverges for lL>n+Q = 
0, i.e. for I = 46.38° or 106.85°; 5p diverges for u N = 0, i.e. for I = sin" 1 (2/ v / 5) = 63.43° 
or 116.57°; and <5p_ diverges for uj n — tl = 0, ie. for 7 = 73.15° or 133.62°. (Note 
that the special value I = sin _1 (2/\/5) coincides with the well-known Newtonian critical 



inclination for J 2 -effects p5| .) As in the equatorial case discussed above, the price to pay 
for taking advantage of the small divisors arising near these poles is the need for a long 
observation period, inversely proportional to the corresponding small divisor. Indeed, the 
dephasing periods between the Newtonian eccentricity vector e^r and the contributions e±,eo 
proportional to ai read 
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ujn 4 — 5 sin I 

where 



r± = — ^ = , (38a) 

u N ±Q 4-5sin 2 /^2cosJ 

T = — = 1 r • 2 r ' ( 38b ) 



B = 



a \ 7/2 /Gm e y 1/2 8tt # e 
a/R^) 7 / 2 x 0.1974 yr. (39) 



[Near each resonance cl>jv ±fi = 0, one of the vectors e± becomes infinitely large with respect 
to the other, and the e + + e_ ellipse degenerates to a circle described with angular velocity 

In order to compare the merits of the different resonances, it is useful to define as a figure 
of merit the ratio of the corresponding perigee displacement Sp by the typical (minimal) time 
T/4 needed for the observation: 

45 p± A l±cos/ 

— — = 4 — cos o 

T± B 2 

R 

= ax — x 1.073 x 10 5 cm.yr" 1 for I = 46.38° and 133.62° , (40a) 
a 

= oti — x 4.509 x 10 4 cm.yr" 1 for I = 73.15° and 106.85° , (40b) 
a 



—— — = 4 — sin o sin 1 
T B 
p 

= a x — x 1.395 x 10 4 cm.yr" 1 for / = 63.43° and 116.57° . (40c) 
a 

These results underline that the inclinations around 46.38° or 133.62° have the capability 
to give the largest effects in a given observational time, as confirmed by the width of the 
resonances in Fig. 2. Equations (|40| ) also show that low orbits (say Rq < a < 2i? e ) seem 
preferable in that they give larger Sp's in a given observational time. [It is easy to see that 
the a -1 dependence of Sp/T still applies for high, tidally perturbed orbits. This justifies 
our concentrating on low orbits, with negligible tidal effects.] However, equations (|3"8"D , (|3"9"D 
show that the inclination / must be fixed with a high precision, typically within less than 
one arc-minute if one wishes to make the fullest use of the observational time T/4 which, 
for practical reasons, will probably not exceed ~ 10 years. The needed inclination can be 
computed thanks to the following asymptotic formulae: 

|/ - /pole | w ' 3.19° ^ for / pole = 46.38° or 133.62° , (41a) 



\ V2 i 

|/-/ P oie|~(^ 2.41°-^ for J pole = 73.15° or 106.85° , (41b) 

\ 7/2 -1 

|/ - /poid « ( ^ 2.83° for 7 pole = 63.43° or 116.57° . (41c) 
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For instance a low satellite (a ~ R®) will give rise to a perigee displacement Sp± ~ «i x 
10 6 cm in T±/4 = 10 yr if the inclination differs from the pole value 43.38° (or 133.62°) by 
only 3.19°/40 = 4'47". An error of one arc-minute on I would change the observational time 
and the perigee displacement by factors ~ (0.8) ±:L . 

The six resonant values of the inclination are solutions of the simple trigonometric equa- 
tions 4 — 5 sin 2 1 + 2x cos 1 = 0, where x G { — 1,0, 1}, and do not depend on any experimental 
data. They will therefore enhance preferred-frame effects in the motion of the natural satel- 
lites of the different planets in the solar system (although tidal forces may not be negligible 
is some cases). However, none of the known natural satellites' orbits has an inclination close 
enough to one of the six poles. On the other hand, there are hundreds of artificial Earth 
satellites, and many of them would allow one to tighten the present bounds (p|),(|ToD on a± 
if the evolution of their perigee could be tracked at the centimeter level. For instance, the 
classes of satellites 'GPS BIT and 'GOES' would typically allow one to measure perigee ra- 
dial displacements of oti x 3 x 10 5 cm in an observational time of T/4 ~ 15 yr. Even if these 
satellites could be tracked with sufficient precision, it would remain to see whether the effect 
of non-gravitational forces would allow one to make full use of such long data span. However, 
a statistical study of these ~ 20 satellites could probably allow a significant reduction of the 
sources of errors. The same remarks can be formulated about the classes of satellites 'GPS' 
and 'Glonass' whose perigees are deviated by ~ ol\ x 1.5 x 10 5 cm in T/4 ~ 8 yr. [There are 
also many low {i.e. fast) satellites, like 'Starlette' or 'Seasat', which give perigee oscillations 
of ~ a% x 3 x 10 4 cm in less than one year, but the air drag is very large for such low 
satellites.] The first Laser Geodynamical Satellite (LAGEOS I) is an interesting candidate, 
not only because its position has been laser tracked for years at the few centimeter level and 
because it is submitted to very small air drag, but also because several contributions to the 
eccentricity vector e are large in its case: e + gives a perigee deviation of «i x 4.2 x 10 4 cm 
in T + /4 ~ 1.9 years, e_ a deviation of a± x 2.0 x 10 4 cm in T_/4 ~ 5.3 months, and eo a 
deviation of a± x 8.8 x 10 3 cm in To/4 ~ 1.2 years. The superposition of these three effects 
thus yields a complex signal which should, hopefully, be distinguishable from other, New- 
tonian, contributions. Unfortunately, the variations of the eccentricity e of this satellite are 



not very precisely modelizable, as explained in Refs. p9|j20|] : it undergoes small oscillations 
of ±3 x 10 -7 , which correspond to perigee deviations of ±4 m. The recently launched second 
LAGEOS, which has an even more favorable inclination, might provide another interesting 
experimental probe of preferred- frame effects. Indeed, a duration of only T + /4 ~ 1.3 yr 
would suffice to observe a radial displacement of its perigee by 5p + ~ a\ x 7 x 10 4 cm. 
The best tool for tightening the present bounds (|9D,([T0|) on a\ would be a (naturally or 
artificially) drag-free satellite, launched on a favorably inclined orbit and tracked for years 
at the centimeter level (via laser ranging or onboard GPS receivers). Let us note that there 
are plans for launching in the near future drag-free satellites with GPS receivers onboard: 
Gravity Probe B (GPB) and the Satellite Test of the Equivalence Principle (STEP). 

For completeness, we now briefly discuss how the above solution (|35|) for the eccentricity 
vector e is modified when tidal forces (which are very small for low satellites) are taken into 
account. The first (trivial) modification is of course that the disturbing function fll5|) now 
gives a contribution to the precession velocity (O) ( |12f| ) of the orbital plane: 
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(O) = -- n cos/ [2J 2 {R e /af + iV 2 /n 2 ] + 0(a ie ) + 0(e 2 ) . (42) 

On the other hand, the effect of tidal forces in equations ( |12b[ ) and ( |12ep for e and uo 
is somewhat more involved [even after averaging over monthly and yearly frequencies as 
explained in section II, Eq. (jl5j)]. The vectorial equation satisfied by the eccentricity vector 
e in the orbital plane reads now 

(d'e/dt) = cx(/ie + kxc) + z/bx(axe) + 0(a x e 2 ) + 0(e 3 ) , (43) 

where 

/x = ~ n [J 2 (R (B /a) 2 (A - 5 sin 2 /) + 2N 2 /n 2 } (44) 

is the analog of the angular frequency un of Eq. (|^), and 

v = — N 2 sin 2 / (45) 

is an additional contribution due to tidal forces. The ratio u/fi can become positive and 
larger than unity when the inclination / is very close to the well-known ^-critical inclination 
/ = sin- 1 (2/ v / 5) « 63.43° or 116.57° (e.g. \I - I \ < 4' for a ~ 2/? ffi , and |/ - / | < 8" for 
a ~ /?©). Equation (|43"D then formally exhibits an exponential blow up of the Newtonian 
eccentricity vector. This indicates that our simplified (linearized and time-averaged) treat- 
ment of the evolution of e becomes inadequate. In the following, we restrict our attention 
to the generic case v f [i < 1. One finds that e^r has no longer a constant norm, as opposed 
to Eq. ( p3|) and (^) above, but that it is now moving on an ellipse, whose axes are directed 
along a and b and have a ratio e a /eb = (1 — v/ fa) 1 / 2 : 

&n = e a cos 1 — v I ii t + est j a + t h sin ^ua/ 1 — u/ ft t + est j b. (46) 

[Note that this ellipse reduces to a circle in the particular case of equatorial orbits (/ = 0) 
that we considered in section II.] The solution of equation ([43D can then be obtained in a 
geometrical manner similar to the one used in Eqs. (|32|)-(p5|). Let denote a generalization 
of the Newtonian eccentricity vector when the frequency [i is replaced by an arbitrary 
constant /ij such that v j \ii < 1, i.e. fii < or /ij > v (we shall see below that these values 
of Hi are sufficient for our purpose). Like e^r in Eq. (|46| ) above, is moving on an ellipse 
in the orbital plane, with a non-uniform velocity which has the sign of /ij, and such that its 
components along a and b are oscillating with a constant frequency ^(l — u/fii) 1 ^ 2 and have 
an amplitude ratio e\j ' e\ = (1 — vjn^l 2 . [The particular cases of constant vectors along 
a or b are obtained respectively for ^ — > 0~ and fa = v]. The time derivative of a linear 
combination e = e N + J2 i ej then reads 

d'e/dt = c x Le| V, (fa ~~ A 4 ) e « + ^ b x (a x e) . (47) 



This has precisely the form of Eq. ([43]) , with k x c having been decomposed as a sum of 
vectors (fa — fa) moving on ellipses as described above. It is easy to check that any vector 
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having a periodic motion in the orbital plane admits such a decomposition, and that the 
corresponding /ij's always satisfy the condition vj\ii < 1. In particular, one finds that the 
source term k x c can be decomposed as in Eq. (|33|), where K is still given by Eq. (|3^ 
but where K ± have more complicated expressions 



K± = cos 5 



+ cos/ 
2 - v/n T 



with 



sin(f2 — a) a + ^ cos(f2 — a) b 



(48) 



/i± = T^V 1 + (^/2^) 2 + v/2 



(49) 



Note that these expressions for K± reduce to those of Eq. (|34a|) when v = 0, i.e. when tidal 



forces are neglected or for equatorial orbits. Hence the solution of Eq. (43) can be written 
simply as 



ejv + «i 



4a 2 c 2 



w 



K 



+ 



+ 



K 



K 



/i /i — /i_ 



(50) 



It is therefore of the same kind as Eq. (p5|) , and there still exist poles for some particular 
values of the inclination, only slightly modified with respect to those of Fig. 2 for relatively 
low satellites. [For the same reason indicated above, our treatment becomes inadequate 
when jj, — > + , formally corresponding to (1 — u/ ji) 1 ^ 2 becoming large and pure imaginary] 



V. PERTURBATIONS DUE TO THE ORBITAL VELOCITY OF THE EARTH 

AROUND THE SUN 

A. Perturbations of the eccentricity vector 

In the previous sections, we have neglected the orbital velocity v ffi of the Earth around 
the Sun, but we shall see below that it can also lead to significant preferred-frame effects 
on artificial satellites, although they will typically be 12 (^ w/v & ) times smaller than those 
proportional to w. The interest of effects involving v e is twofold. First of all, we shall 
see that arbitrarily small divisors can enhance the preferred-frame effects we are studying 
for any value of the inclination / of the satellite's orbit (if the semi-major axis a is chosen 
appropriately) instead of the discrete resonances displayed in Fig. 2. Moreover, the existence 
of a motion around the Sun (with well-defined amplitude and phase) is known for sure 
to be part of the 'absolute' velocities entering preferred-frame effects. By contrast, the 
identification of w with our velocity with respect to the cosmic microwave background is a 
specific assumption. Although this assumption is plausible on field-theoretical grounds, the 
bounds derived on a± in the literature PJl^l and the discussions of the previous sections 
of the present paper are strongly dependent on it. It seems therefore of importance to 
determine limits on a\ by relying only on an unambiguously present velocity such as v e . 

It is straightforward to generalize the results of the previous section to the case of a 
velocity v e which is not constant, as opposed to w. The equation satisfied by the eccentricity 
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vector e is Eq. (pOD, but we concentrate now on the effects generated by the source term 
k® x c, where given by Eq. (p2|), is a vector rotating with angular frequency n® = 
(Gm e /r^ ) 1 '' 2 . As in the previous section, it suffices to decompose k® x c as a sum of 
constant-norm vectors rotating with constant angular frequencies to derive the contributions 
to e involving u®. Let us denote as /® = 23.5° the inclination of the ecliptic with respect 
to the Earth's equatorial plane, and choose the origin of time at the vernal equinox. The 
source term k® x c can then easily be decomposed as 

K ® x c = ""I ^$ v ® ( K ++ + K +- + K + + K ~ + K °+ + K °-) > ( 51 ) 



where 



K ± = — - sin J sin/® [cos(n®t) a=|=sin(n®t) b] , (52) 

and if s, s' denote two independent signs (s = ±1, s' = ±1) 

T ^ / 1 + scosl 1 — ss'cos/® r . . . . , . . . , s 

K ss > = s |cos(si2 + s n®t) a — sm(si2 + s n©r) bj . (53) 

The vectors Ko± are rotating with angular velocities =Fn® in the orbital plane, whereas 
K ss / are rotating with angular velocities — (sQ + s'n©) in this plane. These notations are 
chosen to simplify the expression of the eccentricity vector e written below; it generalizes the 
notations K±,K introduced in section IV, which would be denoted K± ,K o in the present 
convention. The eccentricity vector can then be immediately written as 

Gm ffl / K++ K 

e = e N + at 



4a 2 c 2 \Co N + fl + n® ujN + d — n^ 

K + K__ K 0+ K _ 



ujn — + ™® Wat — — n© + ujn — n f 



(54) 



to which should be added the contributions proportional to w of Eq. (p5|). The solution (|54[) 
now exhibits resonances for any value of the inclination / provided the semi-major axis a is 
chosen appropriately. Figure 3 displays the values of a and / for which one of the divisors 
in Eq. ([54]) vanishes. Note that when a tends formally towards 0, these curves tend towards 
the six poles found in section IV, because u>n and Q [Eqs. ( |3"TD and Q29D] are proportional 
to a _7//2 whereas n e = 2tt/(1 yr) is constant. 

As in the previous section, the price to pay for a small divisor is the need for a cor- 
respondingly long time of observation, say T/4 = ii/2{l)n ± Cl ± n®), or tt/2(lun ± ^®), 
depending upon the concerned divisor. The figure of merit 8p/{T /A) of the different orbits 

can be computed like in Eqs. (0), and Figure 4 shows that the contributions e_j and e h 

give the best results, i.e. the largest perigee displacements in a given observational time. 
This is due to the rather small value of the inclination J® = 23.5° of the ecliptic with respect 

to the equatorial plane, since e H and e_ + are precisely the only contributions which do not 

vanish as J® — >• 0. The values of a and / for which these two main contributions diverge 
have been plotted in solid lines in Fig. 3 in order to distinguish them from those involving 
sin/® or 1 — cos/® (vanishing when /® — > 0) which have been plotted in dashed lines. 
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The largest figure of merit displayed in Fig. 4 is obtained for a m i2® and / = 43.22°, 
and reads 45p + _/T + _ = a\ x 8.62 x 10 3 cm.yr -1 . As expected, it is (w/vq ~) 12 times 
smaller than the best result ( |40a| ) of section IV, which was obtained for I = 46.38° (or 
133.62°). However, the high-precision observation of a low orbit with an inclination close to 
43.22° could lead to the first experimental results about preferred-frame effects free from any 
assumption about the fundamental rest frame. [The satellite LACE happens to be precisely 
on such an orbit, but it is not drag-free.] 

We have seen in section IV that all of the three contributions e + , e_ and e (proportional 
to w) to the eccentricity vector e are a priori large for the orbit of LAGEOS I, thus giving 
rise to a complex signal in the evolution of the eccentricity vector. It is interesting to quote 
that the contribution e_ + (proportional to v®) is also enhanced by a small divisor for this 

satellite, since the amplitude 5p |_ is about oc\ x 2 x 10 3 cm in a typical observational time 

of T_ + /4 ~ 7 months. It would be therefore interesting to analyze the LAGEOS I data to 
look for preferred-frame signals of both the w and the v e types [if it turns out feasible to 
extract such information in presence of the non-modelizable forces that act on the satellite]. 

A remarkable feature of Fig. 3 is the existence of maxima in the curves representing the 
loci of the v ffi -type resonances in the (I, a) plane. These maxima occur at / = or 180° 
(for a/R® = 1.94, 2.36 and 2.65), / = 78.46° or 101.54° (for a/R e = 1.67), and I = 90° (for 
a/R® = 1.59). The amplitude of the v e -preferred-frame effects for orbits located close to 
these values of I and a would be therefore almost insensitive to injection errors or fluctuations 
in I. This suggests that these orbits might be especially robust tools for constraining this 
type of effects. The best choices would be the orbits located close to (J = 0, a = 1.94i? e ) 
or (/ = 101.54°, a = 1.67i?©), which correspond to rather large figures of merit 4Sp/aiT in 
Fig. 4 (respectively 5.15 x 10 3 cm.yr -1 and 3.58 x 10 3 cm.yr -1 ). 

B. Along-track perturbations of the satellite 

We have not yet discussed the perturbation of the element a, which is related to the 
angular position of the satellite at t — 0. Using the expression of the disturbing function 
(|TB|) in ( |12d| ), the secular evolution of a is given by 

/ da\ . ai n 2 a [w + v e ,c,e] am, , 2 . , , , 

\Tt) = aN + T — ? ^ (w + v e ) + 0(a lC ) , (55) 

where &n is the Newtonian contribution due principally to tidal forces (|i~5|) and to the Earth's 
quadrupolar moment (|T4l). The second term of the right-hand-side, though proportional to 
1/e, does not lead to any interesting effect in the limit of a small eccentricity. Indeed, it 
precisely cancels a corresponding term in the secular evolution of to. In other words, the 
angular position of the satellite with respect to the ascending node is not affected by this 
type of contribution. 

By contrast, the last term proportional to (w + v e ) 2 = w 2 + + 2(w ■ v e ) induces 
interesting effects on the position of the satellite through the time variation of w • v ffl . [The 
secular drift in a due to the constant term w 2 + v 2 ^ is unobservable because, as we shall see 
below, it can be absorbed in a small renormalization of Newton's constant G.} Indeed, let 
us denote as before the right ascension and declination of w as (a, 5), the orbital frequency 
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of the Earth as n©, and the inclination of the ecliptic as 7®. Equation (|55|) can then be 
rewritten as 



(a) = est — 2ot\n 



1 + COS la 



cos 5 sin(n®t — a) + 0(1 — cos/®) + 0(sin/ g 



(56) 



where we have written down only the largest contribution. The integration of this equation 
shows that the longitude of the satellite is modulated by an oscillating term 



5a = at 



n W Va 



n® c c 



'1 + cos/ ffi ) cos 5 cos(n ffi t — a) 



corresponding to an along-track oscillation 



5x/i = a5a « a\(R^/ a) 1 ' 2 cos(n®t — a) x 9.17 x 10 5 cm. 



(57) 



(58) 



The large amplitude, yearly periodicity, and weak altitude dependence (oc a -1 / 2 ) of this effect 
makes it a promising way of improving the precision of measurement of oti (maybe down to 
the 10 -6 level). However, detailed numerical simulations are needed to assert whether the 
effect (158|) can be separated from the other yearly perturbations. 

Let us note finally that the results (|5T|),(|5"8"|) can be derived in a totally different way, 
by starting directly from the Lagrangians ([]]) and (|TTD . The terms proportional to ct\ can 
indeed be interpreted as a time-dependent renormalization of the gravitational constant 
experienced by the satellite 
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G{t) = G 



at 
2d 1 



w 



(59) 



Taking into account the adiabatic constancy of the Delaunay variable (or "action variable" ) 
L = [G(t)m©a(t)] 1/ ' 2 , Eq. (|59|) induces a variation of the semi-major axis 



a(t) = a N 



2^ (W + Ve) 



-i 



(60) 



In the limiting case of small eccentricity, the angular position 9 of the satellite with respect 
to the perigee can be identified with the mean anomaly £ — J n(t) dt + aN, where aN is a 
constant and where the orbital frequency n(t) reads 



n(t) 



( G(t)m, 



1/2 



n N 



1 - 



2? 



; w + v, 



(61) 



Integrating Eq. (|6lD reproduces the result (|57|). Note that here again a very small instanta- 
neous perturbation 0(v 2 /c 2 ) has been enhanced by a large factor n/n®. As an aside, let us 
remark that another consequence of the time-dependent renormalization (|59| ) of the gravi- 
tational constant is to cause a yearly "breathing" of the radius of the Earth, with associated 
yearly variations of its moment of inertia and of its angular velocity (see 0). However, 
the amplitude of these variations, given the present limits on a\, are too small to be of 
observational significance, e.g. SR^it) = — (din R^/dliaG) aiR® w • v e (t)/c 2 < 0.4 mm. 
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VI. CONCLUSIONS 



Artificial Earth satellites can be very useful tools to probe the field content of gravity, i.e. 
specifically whether it contains a vector or second tensor interaction leading to preferred- 
frame effects in local gravitating systems. Thanks to the appearance of small divisors which 
enhance the preferred-frame effects on the eccentricity vector when the inclination and/or 
the semi-major axis of an orbit are chosen appropriately, it seems conceivable to tighten 
the present experimental bounds on the preferred-frame parameter a% down to the 10~ 5 — 
10~ 6 level. What is needed are centimeter-level tracking data of a (naturally or artificially) 
drag-free satellite over time scales large enough to separate from Newtonian contributions 
the secular motion of the eccentricity vector induced by a non-zero ct\. Among the class 
of zero-inclination (equatorial) orbits, geostationary ones are nearly optimal. Along-track 
oscillations of a satellite with yearly period constitute another promising way of measuring 
ct\ around the 10 -6 level. There is also a wide class of orbits for which preferred-frame 
effects on the eccentricity vector due to the orbital velocity of the Earth around the Sun 
are enhanced by small divisors. These could be used to obtain the first bounds on a.\ 
independent of any hypothesis concerning the gravitationally preferred rest frame. 

In this paper, we have presented an approximate analytical treatment of the dominant 
preferred-frame effects. The value of this treatment is mainly indicative, as a help for select- 
ing the most favorable orbits. In practice, we advise to resort to direct numerical integration 
of the equations of motion (and, evidently, to a multiparameter fit to the experimental data). 
For the convenience of the interested reader, we end by giving the ol\ (and for completeness 
a<i) contributions to the equations of motion. Contrary to the rest of the paper, we have in 
mind here global (barycentric) equations of motion (written in a post-Newtonian coordinate 
system appropriate to the description of the entire solar system; see e.g. 
acceleration of a satellite with respect to the Earth has the form 
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where the general relativistic geocentric acceleration (including a relativistic treatment of 



multipolar and tidal effects) will be found in full detail in Ref. where A non - grav 
Fnon-grav/^sat denotes the acceleration induced by non-gravitational forces, and where ( 
lowest order in the deviation from general relativity, and for m sat <C m®, r <C r e0 ) 
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In these formulae, r and v denote the radial distance and the velocity of the satellite with 
respect to the Earth, n the unitary vector directed from the Earth to the satellite, n e0 the 
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unitary vector directed from the Earth to the Sun p3| , and w e = w+v e the absolute velocity 



of the Earth with respect to the gravitationally preferred rest frame, v e being its orbital 
velocity around the Sun. For conceptual clarity we have indicated within square brackets 
the contributions of a% and a 2 due to the violation of the strong equivalence principle, i.e. 
the terms generated by expanding the //-dependent term coming from the Lagrangian Lp^ ir] , 
Eq. (§). The different terms of Eqs. (|63|) are classified by order of decreasing magnitude. 
[The even smaller contributions due to the coupling of the satellite to the Earth's intrinsic 
angular momentum are given in Eq. (9.20) of Ref. ||; they cause an additional secular 
precession of the satellite's orbit which is independent of the absolute velocity with respect 
to the preferred frame.] 

The contribution proportional to u> e n in Eq. fl63a|) is responsible for the along-track 



perturbations of the satellite studied in section V.B above, whereas the one proportional to 
(nx w e ) x v, together with the contribution proportional to (Gm^/r) n eQ in the second line, 
is responsible for the perigee displacements studied in the rest of the paper. Note that the 
(Gm^/r) n eQ term has the same form as the one due to a violation of the strong equivalence 
principle [with an opposite sign: 2E| rav /m e c 2 « —9.2 x 10~ 10 , Gm^/ac 2 « (i?©/a) x 7.0 x 
10~ 10 ]. Actually, neither of these two terms is of much observational significance, as their 
contributions are much smaller than that due to the (n x w e ) x v term. 

The force proportional to (n e0 x w) x v in the second line of Eq. ( |63a|) has not been 
considered in the rest of the paper. It arises in the connection between the locally inertial 
geocentric frame used in the body of the paper and the global barycentric one used here, 
and has the form of a Coriolis force. This Coriolis force acts also on gyroscopes (including 
the spinning Earth) and adds up to several other relativistic effects causing a universal 
precession O of gyroscopes and satellite's orbits with respect to the barycentric frame (see 
section 9.1 of |J and [0). The «i-dependent contribution to O reads 

Gm & . s 

°«i = ai n ®0 x w • ( 64 ) 

Although the time dependence of Q ai , via n ffi0 (t), might in principle allow one to separate 
it from the other relativistic contributions to f2 (when discussing observables related to the 
global barycentric frame), its magnitude is too small to be of observational significance. 
[It yields a yearly oscillation of the satellite with respect to the barycentric frame which is 
smaller than (a/i? e ) x 0.1 mm.] 
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FIGURES 



FIG. 1. Definition of the orbital elements e, I, f2, oj, and of the rotating orthonormal basis 
(a, b = c x a, c) linked with the orbital plane. In the limit of a small eccentricity e, the angle 
between the perigee and the position of the satellite at t = can be identified with the orbital 
element a. 

FIG. 2. Expected radial displacements of the perigee as functions of the inclination of the 
orbital plane. The solid line correspond to the contribution 6p + , the dashed line to Sp-, and the 
dotted line to 5po- 

FIG. 3. Values of the semi-major axis of the satellite's orbit for which the displacement 5p of 
the perigee has a resonance, as functions of the inclination of the orbital plane. The bold lines 
correspond to the six poles shown on Fig. 2, which exist for any value of a (small enough for 
tidal forces to be negligible). The solid lines correspond to effects proportional to which do 
not vanish when the inclination i® of the ecliptic with respect to the equatorial plane is neglected. 
The dashed lines correspond to the vanishing effects when J® — > 0. The dotted lines correspond to 
unphysical values of a < R® . 

FIG. 4. Figure of merit of the resonant preferred-frame effects corresponding to the orbits of 
Fig. 3, i.e. ratio of the expected radial displacement Sp of the perigee by the typical observational 
time T/4. The plain lines correspond to the contribution 5p+±, the dashed lines to <5p_±, and the 
dotted lines to 5po±. 
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